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Abstract. In this paper we will show how the results found in [7], about the Busemann functions 
in last-passage percolation, can be used to calculate the asymptotic distribution of the speed of 
a single second class particle starting from an arbitrary deterministic configuration which has a 
rarefaction fan, in either the totally asymetric exclusion process, or the Hammersley interacting 
particle process. The method will be to use the well known last-passage percolation description of 
the exclusion process and of the Hammersley process, and then the well known connection between 
second class particles and competition interfaces. 



1. Introduction 

The macroscopic behavior of the density profile of the totally asymmetric nearest neighbor ex- 
clusion process (TASEP) is governed by the Burgers equation, which corresponds to the "shape 
theorem" in last-passage percolation [17 . The second class particles, that follow roughly the be- 
havior of a perturbation of the system, are asymptotically governed by the characteristics of the 
Burgers equation. When there is only one characteristic, the second class particle follows it; when 
there are infinitely many (rarefaction fan), the particle chooses one of them at random to follow. 
Until recently, there was only one method, introduced by Kipnis and Ferrari [10J, to compute the 
distribution of the asymptotic speed when the initial distribution is a product measure with den- 
sities p £ [0, 1) and p' G (p, 1], to the right and left of the origin, respectively. They proved that 
the rescaled position of the second class particle converges in distribution, as time goes to infinity, 
to a random variable uniformly distributed in the interval [1 — 2p' , 1 — 2p\. Their method is based 
on the hydrodynamic description of the macroscopic behavior and it strongly relies on the product 
structure of the initial profile. 

Later, following initial results by Guiol and Mountford [13J and Ferrari and Pimentel [12j . Ferrari, 
Martin and Pimentel in [11] proved that for initial configurations exhibiting a rarefaction fan, the 
second class particle almost surely have an asymptotic speed, but no new results were found about 
the distribution of this speed. 

Recently, in a paper by Amir, Angel and Valko [2], a new method has been introduced to 
study the behavior of a multiclass TASEP, which allows for the calculation of joint distributions of 
multiclass particles under a specific initial configuration. However, as far as we can see, this does 
not shed new light on how to calculate the asymptotic distribution of a single second class particle 
in an arbitrary initial configuration exhibiting a rarefaction fan. 

We will also consider the Hammersley interacting particle process [I], where the situation con- 
cerning the second class particle in a rarefaction fan is very similar to TASEP: almost sure existence 
of the asymptotic speed is proved by Coletti and Pimentel in [9] (they even prove this for more 
general objects than second class particles), but the distribution of the second class particle in a 
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rarefaction fan was only determined for a very specific family of initial conditions, involving Poisson 
processes ([HE]). 

In this paper we will show how the results found in [7J about the Busemann functions in last- 
passage percolation (LPP) can be used to calculate the asymptotic law of the speed of a single second 
class particle starting from an arbitrary deterministic configuration which has a rarefaction fan, in 
either the totally asymetric exclusion process (TASEP), or the Hammersley process. The method 
will be to use the well known LPP description of the TASEP and of the Hammersley process, and 
then the well known connection between second class particles and competition interfaces (see |12| 
for the TASEP case, and [7J for the Hammersley case) to show that the limit law can be expressed 
in terms of the supremum of two independent random walks in an environment that depends upon 
the initial configuration. This method allows us to recover the classical results, as well as to get 
explicit formulas for the law of the asymptotic speed starting from periodic profiles. Busemann 
functions in First Passage Percolation were first introduced by Newman [14J, and these were later 
extended to LPP by Wiithrich in [18J, by Ferrari and Pimentel in [12j and by the authors in [7J. 

To give some intuition for the main idea in the paper, we present a heuristic outline of it in the 
lattice LPP competition model. Suppose that to each x G Z 2 we attach an exponential random 
variable X x of intensity one. For x < y (coordinate- wise) let L (x, y) denote the maximum, over 
all up-right nearest-neighbor paths w connecting x to y, of the sum of X z along zu, leaving out x 
(see (I2.6p ). We adopt the notation in [7J and, to indicate that it refers to the lattice model, we use 
the superscript £. Now choose two non-ordered points x = {x\,X2} and y = (j/i , 2/2) with y\ < x\ 
and i)2 > X2- We will assign each point z S R 2 with x,y < z to either x or y in the following 
way: x "conquers" z if L (x, z) > L (y,z), and otherwise y conquers z. The point z gets the 
color blue if it is conquered by x and otherwise it gets the color red. This example corresponds 
to the two-corner profile discussed in Section 14.11 The competition interface is the up right path 
that separates the blue and the red regions. It has been shown in [12] that a competition interface 
has an (a.s.) asymptotic angle 0, which will be random in this case (rarefaction regime). To check 
whether < a, for some fixed a G (0, vr/2), one can consider z Q far away on the line with angle a, 
and compare the length L (x, z a ) with l/(y, z Q ). It is at this point that the Busemann function 
comes into play: for a G (0,7r/2), 

a.s. 3 B e a (x,y) := lim L e (y, z a ) - L e (x, z «) • 

This means that checking whether i/(x, z a ) < l/(y, z Q ) (this corresponds to < a) eventually, 
corresponds to checking whether I?^(x,y) > Q- In [7j we derived the distribution of B^, which 
is why we are in fact able to calculate the asymptotic distribution of 0. This argument can be 
extended for any deterministic initial profile in the rarefaction regime. 

Overview. In Section [2] we state the main results, Theorem [U Theorem [2] and Theorem [3] . In 
Section [3] we introduce the equilibrium description of the Busemann functions and establish its con- 
nection with second class particles and competition interfaces to prove all the theorems. Section [4] 
and Section [5] are devoted to compute the law of the asymptotic speed, for some initial configura- 
tions in the rarefaction regime, in the TASEP and in the Hammersley models, respectively. 

2. Main Results 

In this section we state the general results about the asymptotic distribution of the second class 
particle in a rarefaction fan in the TASEP case and the Hammersley case. The TASEP result is 
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formulated in two ways: the first using the particle description, the second using the last passage 
description. In Section [3] we will link the two results. Applications of these general results to 
particular initial conditions can be found in Section [J] and Section [5j 

2.1. Second-class particles in the exclusion process. The one-dimensional nearest neighbor 
totally asymmetric simple exclusion process is a Markov process (r^ , t > 0) in the state space 
{0, 1} Z . In our notation, rj t (x) = 1 indicates that there is a particle at site x at time t; only one 
particle is allowed at each site. At rate one, if there is a particle at site x € Z, it attempts to 
jump to x + 1; if there is no particle in x + 1 the jump occurs, otherwise nothing happens. To 
construct a realization of this process, one considers independent one dimensional Poisson processes 
M = (N x (-),x G Z) of intensity 1. The process (rj t , t > 0) can be constructed as a deterministic 
function of the initial configuration 77 and the Poisson processes M as follows: if s is a Poisson 
epoch of N x and there is a particle at x and no particle at x + 1 in the configuration t/ s _ , then 
at time s the new configuration is obtained by making the particle jump from x to x + 1. Let $ 
be the function that takes 77 and A" to (r/t , t > 0). Let 77 and 7/ be two arbitrary configurations. 
The frasic coupling between two exclusion processes with initial configurations 77 and 77' respectively 
is the joint realization ($(77, A/"), 3>(f/> A/ - )) = ((r]t,r]' t ), t > 0) obtained by using the same Poisson 
epochs for the two different initial conditions. Liggett (1985, 1999) are the default references for 
the exclusion process. 

Let 77 be an initial configuration and let 77' be the configuration that differs from 77 only at site 0. 
With the basic coupling, the configurations at time t > differ only at one site A r '(t) defined by 

X\t) :=^xl{77 t (x)^(z)} 

(1 denotes the indicator function). The process (A' ? (t), t > 0) is the trajectory of a "second class 
particle". This particle jumps one unit to the right at rate one if there is no 77 particle in its right 
nearest neighbor and it jumps one unit to the left at rate one if there is an rj particle in its left 
nearest neighbor site, interchanging positions with it. 

To state our results, we need some more notation. For each 77 we associate a new configuration 
77 as follows: 



(2.1) fj(k) :- 



rj(k) for k < 
for k = 



1 for k = 1 

rj(k — 1) for k > 1. 

The hole-particle pair at the origin will be important to formulate the LPP competition model and 
to establish the connection between second class particles and competition interfaces [12]. Note 
that 77(0) is not relevant for 77; this corresponds to the idea of having a second class particle in 0. 
For p £ (0,1) let {Exp k (p) : k £ Z} and {Exp' fc (l — p) : k G Z} be independent collections of 
i.i.d. exponential random variables of intensity p and 1 — p, respectively. These collections are also 
assumed to be independent of 77, whenever 77 is random. For k > let 

(0 9 x yP ,v ._ / - Ex Pfc(p) if V( k ) = 

[ ' k 1 Exp' fc (l-P) if77(fc) = l, 

and for k < let 

( o\ yP ,ri ._ / Ex Pfc (p) if fj(k) = 

[ ' k I -ExpUl-p) if77(A0 = l. 
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Define the random walks 

n —n 

S^ + :=^Zr^dS^:=^Zr. 

k=l k=0 

We assume that the particle system is in the rarefaction regime, i.e.: 

-l 



1 1 
(2.4) p„ := limsup— > n(k) < liminf — > n(k)=:pL. 

k=l k=—n 

Theorem 1. For each u € [1 — 2p^, 1 — 2p v ] let 



l + u 

Pu = n 



If (E3) holds then 

X"(i) v 



U v as t — > oo . 



where V 1 is a random variable with support [1 — 2p' 1 — 2p v ] and for u £ (1 — 2p' , 1 — 2p r] ) 



(2.5) P (£/" > u) = P sup S^tL" > sup fi^f 

Vn>0 ' n>l 

To give some intuition for this result, note that the asymptotic drift of the random walk S^f 1 is 
lesser than or equal to — p'j (1 — p u ) + (1 — p'^ )/p u , which is strictly negative if u > 1 — 2p' v and it 
is a decreasing function of u. Also, the asymptotic drift of the random walk S^ U J^ is lesser than or 
equal to — (1 — p v ) / ' p u + p v / (1 — p u ), which is strictly negative if u < 1 — 2p r) and it is an increasing 
function of u. Note that when u = 1 — 2p' rj or u = 1 — 2p^, (12. 5p may not hold in general! An 
example of this can be found in Section [3] (see Remark 13. 1 H . The left-hand side of (12.5P is still 
determined at the boundary by left-continuity. 

2.2. Competition interfaces in the lattice last-passage percolation model. Let X := 

{X z } xgZ 2 be a collection of i.i.d. random variables with an exponential distribution of parameter 
one. For x, y € Z 2 , with x < y (coordinate-wise), let 

(2.6) L*(x,y) := max Vl z , 

where w : x y denotes an increasing (or up-right) path connecting x to y, not including x. The 
random variable Z/ is called the last-passage (percolation) time from x to y. 

Using the well known connection between the totally asymmetric exclusion process (TASEP) and 
the last-passage percolation model, introduced by Rost [IT], we identify the initial configuration fj 
of particles and holes, that includes a hole at and a particle at 1, with an initial (growth) profile a 
in Z 2 , in the following way: we define the profile a = (as a function of r/), which is a down-right 
bi-infinite sequence of points in Z 2 containing the origin, such that 

' (1,0) if = 0, 



<t(0) = (0,0) and a(k) - a(k - 1) 



(0,-1) if*7(fc) = l. 



Given such an initial profile a, we can define a competition interface, following [12j and [llj . The 
competition for territory goes as follows: choose t G M 2 such that {x : x < t} n a ^ 0. Define 

L l (a~,t) := sup{L e (a(k),t) : k < and a(k) < t} 

and 

L e (a + ,t) := sup{L e (a(k),t) : k > and a(k) < t}. 
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If {x : x < t} n = 0, we define Lf(a^- , t) = 0. The point t gets the color blue if it is further 
(the longest path is longer) from the negative part of a than from the positive part, so if 

L e {a-,t) > L i {a + ,t), 

and otherwise it gets the color red. Note that because of the hole-particle pair at the origin, <r(0) 
can never be the most "distant" point to t. Also, the coloring is well defined for every t above a 
if the initial condition n has infinitely many particles to the left of the origin, and infinitely many 
holes to the right; we will always assume this. The competition interface ffl := (</>n)n>o is the up 
right path through the lattice points starting at that separates the blue and the red region (see 
Section 2 of jllj for an extensive description). 
Under (I2.4p . the rarefaction condition, let 

arctan | j — \ ) < arctan [ | „ v , | j =: a' . 



Pv 




Note that p v /(l — p v ) is the (negative) asymptotic slope of a at +oo, whereas p' v /(l — p' v ) is the 
slope of a at — oo. 

Theorem 2. For each a G [a^,a^] let 



ret • / • 

1 + ytana 

If (|2.4p holds then the competition interface (</>n)n>o satisfies 
(2.7) — ifrr — > (cos @ v , sin @ v ) as n — > oo , 

\4>n\ 

where Q v is a random variable with support [a^a'] and for a G (a v ,a' 



(2.8) P (G* 7 < a) = P (sup S**!? > sup S p n a ? 

\n>0 ' n>l 



Theorem [2] and Theorem [TJ although each of separate interest, are of course intimately linked 
through the connection between TASEP and the lattice LPP, which enables us to directly translate 
the competition interface in the LPP into the position of the second class particle in TASEP; see 
Lemma [3.5l in Section [3] for the details. Our proof will use the LPP description, as explained in the 
Introduction. Again we have that at the boundary, (|2.8p may fail. 



Remark 2.1. Ferrari, Martin and Pimentel in [TT] show that the competition interface in the 
rarefaction fan has in fact almost surely an asymptotic direction (which implies that the second 
class particle in TASEP has an asymptotic speed). They assume an asymptotic inclination of the 
initial profile a in the positive and in the negative direction (which corresponds to an asymptotic 
particle density left and right of zero in TASEP). However, it is not hard to adapt their proof only 
slightly to get the same results under our rarefaction fan condition (|2.4p . 

2.3. Second class particles in the Hammersley particle system. Let P C I 2 be a two- 
dimensional Poisson random set of intensity one. The last-passage time (or longest increasing 
path) between z, v G M?, with z < v, is defined by 

(2.9) L(z,v) := max \w\ , 

where © : z / v denotes an increasing path from z to v, consisting of points in P that are not 
on the same horizontal nor vertical level as z. Furthermore, \w\ denotes the size (the number of 
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Poisson points) of such a path. The Hammersley model can be seen as a Markov process on the 
space J\f consisting of the set of all counting processes v on M such that 

(2.10) liminf^>0, 

y-s>-oo y 



where 

(2.11) u(y) 



KM) fory>0 
-v((y,0]) fory<0. 



(We allow v(y) = — oo for all y < 0.) For each v G J\f, the process 

(2.12) L u (x,t) := sup {v(z)+L((z,0),(x,t))} for x G M , and i > 

is well defined and the time evolution defined by 

(2.13) M?((x,y]) := L„(y,t) - L„(x,t) , 

is a Markov process in J\f, called the Hammersley interacting particle system pp. One can think of 
Mi as a particle system (or process in M) that starts with the atoms (or particles) of v, and that 
moves according to the (informal) rule: for (x,t) G P, at time t > the nearest particle to the 
right of x is moved to x, with a new particle created at x if no such particle exists. 

To introduce a second class particle in the system, we put an atom at in the process v by 
defining 

v([0,x]) =v([0,x]) + 1 for x>0. 

With this new process, and using the same P, we construct M t ^, and then define the location of 
the second class particle Y"(t) as 

Y v (t) = inf{x > : Mj?(x) = M?(x) + 1} . 

We note that Y"(t) is a non-decreasing function of i, meaning that the second class particle moves 
to the right. We use the notation of [7], where relevant details can also be found, including the 
following important connection between the longest path description and the second class particle. 
Let v + , v~ be the processes defined by 

, ni ., -L- / \ f v(x) for x > , _ , N f for x > 

(2. 14) v + (x) = { K ' , - and v (x) = < , , , ~ _ 

I —oo for x < U. [ v{x) tor x < U. 

Then 

(2.15) {Y u (t) <x} = {L u+ (x, t) > L v - (x, t)} . 

The above equation resembles the situation where one should check if the competition interface 
is to the left or to the right of some angle a. This indicates that the speed of a second class particle 
can be related to the respective Busemann functions in a similar way. We require the following 
rarefaction assumption (recall (|2.1ip and (|2.14p ) 

( x} is (v^) 

(2.16) a^—limsup — ^<liminf y -^=:b v . 

x^oo X y^-co y 

For p > let v p be a one-dimensional Poisson counting process of intensity p > 0, that is indepen- 
dent of v } whenever u is random. 



V u as t -> oo . 
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Theorem 3. For each v G [b~ 2 ,a~ 2 ] let 

1 

If d2JBD ZioWs t/ien 

t 

where V u is a random variable with support [b„ 2 , a~ 2 ] and for v G (b~ 2 ,a~ 2 \ 
(2.17) P(V" < v) = P (supM*) - v pv {z)} > sap{v(z) - v Pv (z)}) . 

\z>0 z<0 J 

Condition (j'2.16p ensures that with probability 1, at least one of the two suprema is finite, and 
hence the right hand side probability is well defined for every v G (b~ 2 ,a~ 2 ). Again we have that 
at the boundary, (|2.17|) may fail. 

3. Proof of the Theorems 

We will start with the proof for the Hammersley process, since the ideas are most clearly repre- 
sented in that case. In the TASEP case, we will use very similar methods. 

3.1. Busemann functions and equilibrium measures in the Hammersley model. In the 

Hammersley last-passage percolation model, Wiithrich [18] proved that, for a G (0, 7r/2), 

a.s. 3B a (x,y):= lim L(y, z a ) - L(x, z a ) . 

|z a |-s>+oo 

Here, z a follows the ray through (cos a, sin a). By the invariance under the map a — > a + tt, it 
is straightforward to see that the analogous result holds for a G (ir,3ir/2). The most important 
property of the Busemann functions is the following connection with the associated interacting 
particle process (Theorem 5.3 and Corollary 5.5 of [7]): Let a G (tt,3tt/2) and define 

u a ((x,y]):=B a {0,(y,0))-B a (0,{x,0))=B a ({x,0),(y,0)) for x<y. 

Then v a is the unique spatially ergodic counting process with intensity 

Ei/ a (l) = an a , 

that satisfies (time invariance) 

(3.1) Mt a = u a , for all t > 0. 

A consequence of this result is that we can compute the law of the Busemann function along the 
horizontal axis. For p > let v p be an one-dimensional Poisson counting process of intensity p > 0. 
Then 

( 3 - 2 ) "a = ^v^JT^ • 

We also need the following lemma for the Hammersley process L u (see (j2.12|) ). This result, in a 
completely different formulation, would also follow from Lemma 2.2 in [5]. 

Lemma 3.1. Consider the Hammersley process L u , and suppose that for t >0,0<x<y< Y u (t). 
Then 

L v (y, t) - L v (x, t) < L(0, (y, t)) - L(0, (x, t)). 
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Proof : The fact that x < y < Y u (t), implies that (see (|2.15p ) there exists t x < t y < such that 

L v (x, t) = L((t x , 0), (x, t)) + v(t x ) and L„(y, t) = L((t y , 0), (y, t)) + v{t y ). 

Define z as the (an) intersection point of the longest upright path from (t y ,0) to (y,t) and the 
longest upright path from to (x,t). Then we can see that 

L u ( X ,t) > L((ty,0),(x,t))+U(ty) 

> L((ty,0),z) + L(z,(x,t)) + v{t y ). 

This implies 

L u (y,t) - L u (x,t) < L((t y ,0),(y,t))-L((t y ,0),z)-L(z,(x,t)) 

= L(z,(y,t))-L(z,(x,t)) 

= L(z,(y,t))-L(0,(x,t)) + L(0,z) 

< L(0,(y,t))-L(0,(x,t)). 

□ 



3.2. Proof of Theorem [3J For each v,t > let a = a v = arctan(l/u), x = xt := t/tana, 
Xq := (i,a;tana). By (|2.15p . 

{Y"(t)<x} = {L u+ (x a ) >L,-(x a )} 

= { sup [v + (z)+L((z,0),x a )} >Bup[v-(z)+L((0,z),x a )] } 

^ 0<z<x z<0 > 

= { sup [u+(z) + L((z,0),x a ) -L(0,x a )] 

L 0<z<x 

(3.3) > sup [v-(z) + L((0, z), x a ) - L(0, x Q )] ) . 

2<0 J 

For any compact set KcR 2 

(3.4) (L(z,x a ) - L(0,x Q ) ; z G if) (B a (0,z) ; z G if) as re -> oo . 

This remark is the core of the proof. It follows from Theorem 2.2 in [TJ: we can take n > 1 big 
enough, so that K C [— n, n] x [— n,n]. The a-rays starting at (—n,n) and (n,—n) will coalesce 
at some point c G M 2 . Furthermore, the longest paths to x Q , starting at (—n,n) and (n, — n), will 
converge to the respective a-rays in a bigger bounded square, containing c. From that time on, for 
all z G K , we will have 

L(z,x a ) - L(0,x Q ) = B a (0,z). 
What follows are some technicalities to show that the location of the maximum behaves well. 

Lemma 3.2. If \J tan a G {a v ,b v ) (recall (I2.16P ) then, as x — ^ oo, 

sup [v + {z) + L((z,0),x a ) - L(0,x a )l sup [^ + (z) - v+{z)] , 

0<Z<X 2>0 

and 

sup + L((z,0),x a ) - L(0,x Q )] sup [^"(z) - . 

z<0 z<0 
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Proof : We will prove that a u < Vtana implies the first statement. The proof of the second 
statement follows exactly the same reasoning. 

Pick a < y / tan(a) and C > such that v + {x) < ax + C for all x > 0. For fixed x we know that 

(3.5) (L(0,x a ) - L((z,0),x Q ) ; < z < x) = (L(0,x a ) - L(0,x a - {z, 0)) ; < z < x) . 

Choose a' and a' such that a < a' = V tan a' < Vtan a. Start a stationary Hammersley system 
with a Poisson process of intensity a', z7 a /, using the same Poisson process in the upper half-plane 
as L. Define as before 

Lp a/ (.x,t) = sup{L((z,0), (x,t)) + i> a '(z)} ■ 

We know that the characteristic direction of this process is a', and from [6] we conclude that 
the second class particle of this process will be to the right of x a with probability greater than 
1 — 0(x~ 3 ); call this event E x . Fix < z < x. Because of Lemma [3. 11 under the event E x , we have 

L(0, x Q ) - L(0, x Q - {z, 0)) > L Pa , (x Q ) - L Pa , (x Q - (z, 0)). 

As a process in z, the right hand side is just a Poisson process with intensity a' . This process will 
be above the line az + C eventually, and therefore be larger than v + (z). Using (I3.5D . this means 
that for each e > 0, we can find R > such that for all x large enough, 

PI sup [z/+(s) + L((z,0),x a )-L(0,x a )] > sup + L((z,0),x Q )-L(0,x Q )] ) < e. 

Y0<^<a: 0<^<K y 

Now we can use fl3.4f) to finish the proof of the first statement of Lemma 13.21 □ 



Notice that, since v a is a function of the underlying two dimensional Poisson random set P, 
it is independent of v. Together with (j3.2|) . (|3.3p and Lemma 13.21 this proves the theorem for 
tana G (a v ,b v ). Now, if b v < Vtana then a.s. sup J/<0 v~(y) — v~(y) = oo. By the same reasoning 
as in the proof of the last lemma, 

sup [V - ^) + L((z, 0),x a ) — L(0,x a )] —> oo , as x — > oo . 

z<0 

Hence, both functions in f|2. are zero. The case where \/tana < follows from a similar 
argument. 

Remark 3.1. Note that the limiting distribution of t~ 1 Y' y (t), for a deterministic v, may only have 
an atom at a such that Vtana = or Vtana = 6^. If this is the case, then the function (|2,17p is 
not continuous at these choices of a, but it is continuous everywhere else. An example would be if 
v+ {y) = \y\ — Wy\ 2 ^\- I n this case, when we take a = 7r/4, which corresponds to a u = Vtan a, 
and we make sure that a u < b u , we see that 

P sup [u + (y) - v a (yj\ > sup [u~(y) - u a (y)] < 1, 
\y>o y<o J 

whereas P(lim t ^. 00 Y v (t)/t < 1/tana) = 1. This means that P(lim t ^ 00 Y v (t)/t = 1) > 0. 
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3.3. Busemann functions and equilibrium measures in the TASEP-LPP context. For 

the lattice last-passage percolation model with exponential weights, the existence of the Busemann 
functions was proved in [12] : for a G (0,7r/2), 

a.s. 3 B e a (x,y) := lim l/(y, z a ) - L e (x, z a ) . 

By the invariance under the map a — > a + tt, it is straightforward to see that the analog result 
holds for a G (tt,3it/2). In this set up one can also define last-passage times L v t similar to (|2,12p . 
with respect to an initial configuration ir defined on Z, and an evolution M similar to (|2.13p . This 
evolution will be a Markov process and 

ui((x,y}) := B e a (0, (y,0))- B e a (0, (x,0)) for x < y , 

will be the unique spatially ergodic process with intensity 

Eui(l) = —^=(= Pa ), 
1 + ytana 

that is time invariant (Theorem 8.1 of [7]). This allow us to compute the distribution of B e as 
follows: 

Lemma 3.3. The Busemann function is additive, 

Bi(x,y) = Bi(x,z) + B£(z ) y) J 

and anti- symmetric, 

Bi(x,y) = -Bi(y,x). 
Suppose a G (0,7r/2). Define = B f a (o~(k — l),a(k)). Then all Z^'s are independent and 



(3.6) Z k 



v / -Exp{p a } ifrj(k) = 
Exp{l-p Q } iffj(k) = l, 



where Expjp} denotes a exponential random variable of intensity p. 

Proof : The first two statements are general properties of the Busemann function and follow from 
Proposition 3.1 of [7J. Property (|3.6p follows from Theorem 8.1 of [7] (see also [3]) together with 
anti-symmetry. □ 



3.4. Proof of Theorem O For each t > let ^(t) be the intersection of the competition interface 
& v = {4>n)n>i with the line {(x,t) : x > 0}. For each a G (0, 7r/2) let x = xt '■= t/ taxi a, 
x a := (x,xtana). Then 

W v (t)<x} = {L e (a + ,x a )>L e (o-,x a )} 

= I sup L e (a(k), x a ) > sup L e (a(k), x a ) > 

fc>0,o-(fc)<x a fc<0,o-(fc)<x Q - 1 



fc>0,o-(fc)<x Q 



L\o-(k), Xa )-L e (0,x c 



L £ (a(fc),x Q )-L £ (0,x Q )] }. 



> sup 

fc<0,tr(fc)<x a 

We now need to control the location of the maximum on the profile. For this we use the following 
lemma: 
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Lemma 3.4. Suppose a < aL. Almost surely, there exists M > such that for all x > 0, 

sup l/(cr(fc), x a ) = sup L e (a(k),x a ). 

k<0 -M<k<0 

Similarly, suppose a > a v . Almost surely, there exists M > such that for all x > 0, 

sup L e (a(k), x Q ) = sup L e (a(k),x a ). 

k>0 0<k<M 

Proof : This result follows more or less directly from Proposition 3.1 of [11]. That result is stated 
for profiles that have an asymptotic direction, whereas we only have (I2.4p . However, the proof in 
the more general case is exactly the same. The proposition states (translated to our notation), that 
if q and a' are chosen such that tana > tana' > tana^, and {zj} is a sequence of points going to 
infinity with asymptotic direction a', then there are only finitely many z^'s that belong to a longest 
path from er + to any of the x Q 's. If we take M + > such that the x-coordinates of those finitely 
many Zj's are smaller than the x-coordinate of o~(M + ), then it is clear that all the longest paths 
from <t+ to any of the x Q 's must start before a(M + ). We apply Proposition 3.1 of [11] again to the 
negative part of a, obtaining M_ > 0. Now put M = max{ M_,M + }. □ 

The proof now follows as in the Hammersley case, using the convergence to the Busemann function 
on a compact set. Therefore, 



(3.7) P(G ?? < a) = P sup B£(0, a{k)) > sup fl£(0, a(k)) 

\k<0 k>0 

and Theorem [2] is a consequence of Lemma 13.31 together with (j3.7|) . for a G (a v ,a'). 

We need an argument for a < a v (the argument for a > a'^ will be analogous). When considering 
L(cr_,x Q ), Lemma 13.41 tells us that almost surely, there exists M > such that 



Choose < p < p v such that 



L(a_,x a ) <L((-M,0),x a ) 
1 

> 1 - p. 



1 + y/tsu a 

This is possible, since a < a v implies that p a > 1 — p v . Define a sequence & + {k) = (k, [— X p k\), 
with A p = p/(l — p). Because of the rarefaction assumption (|2.4p . we know that there are infinitely 
many points on <r+ below a + . We first show that with probability 1, there exists ko > 1 such that 
for k > ko, 

(3.8) B a ((-M,0),or+(fc))>0. 

This follows from the fact that 

B a ((-M,0),(k,[-X p k\)) = B a ((-M,0),(k,0))+B a ((k,0),(k,[-X p k\)). 

From Lemma [3.31 we know that B a ((—M, 0), (k, 0)) is distributed as the sum of k + M independent 
exponentials with expectation —l/p a , whereas B a ((k,0), (k, [— X p k\) is distributed as the sum of 
[Apfc] independent exponentials with expectation 1/(1 — p a ). Now (|3,8p follows for k big enough 
from the observation that 

1 Xn 1 Arj 

— + > -- — + = o. 

Pa 1 - Pa I-P P 

Now choose z G 0"+ such that z < cr+(fe), for some k > ko. Then B a ({— M, 0), z) > 0. Therefore, 
for x large enough, we will have that L(z,x a ) > L((— M, 0), x Q ) > L(cr_,x a ). Since z G <r + , this 
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means that with probability 1, the angle of the competition interface will eventually be greater 
than a. 

3.5. Proof of Theorem [T]. Together with Lemma [331 below. Theorem [2l implies Theorem [Tl This 
lemma introduces the hole-particle representation of the second class particle to make the link with 
the competition interface. 



Lemma 3.5. For u £ (—1, 1) let 

a u := arctan 

Then 



1 \ 2 
1 — u 



1 + u 



F{U V > u) =F(@ V < a u ) . 

Proof : We label the particles sequentially from right to left and the holes from left to right, with 
the convention that the particle at site one and the hole at site zero are both labeled 0. Let Pj(0) 
and Hj(0), j G Z be the positions of the particles and holes respectively at time of the initial 
configuration fj. The position at time t of the jth particle Pj(t) and the itb. hole Hi(t) are functions 
of the variables L(o~, z): at time L(o~, the jth. particle and the ith hole interchange positions. 

Disregarding labels and defining r) t {Pj(t)) = 1, fj t {Hj{t)) = 0, j G Z, the process fjt indeed realizes 
the exclusion dynamics (see |17|). Now, let tq = and 

T n := max{L((j _ ,0 n ),L((j + ,0 n )} . 

Define the process {(I a (t), J a (t)),t> 0) by 

{F(t),r{t)) = <l> n forte [r n ,r n+ i). 

This process is the hole-particle representation of the second class particle. By Proposition 3 of 

{x^t) , t > o) = (r (t) - r{t) , t > o) . 

By using a similar argument as in the proof of Proposition 5 of [12] , one has that 

— K - J —^ — ->■ /(G) as t -)■ oo 

where 

1 — \/tan a 

/(«) := : ~ n ; . 

1 + ytan a 

which implies Lemma 13.51 □ 



4. Computing the law of the asymptotic speed in the TASEP 
4.1. Two-corner initial profile. Choose integers x,y > 1 and let r] = r\ x ^ y be defined as follows: 

if - (x - 1) < k < or k > y 



r lx,y{}^) :- 



1 if < k < y — 1 or < —x. 



For this initial configuration the supremum is attained at deterministic points and it is easy to see 
that 

supS^'l = T X)P and sup = T y ^ p , 

n>0 ' n>0 
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where T X)P and ^ y ,i- P are two independent gamma random variables of parameters (x, p) and 
(y, 1 — p), respectively. On the other hand, 

p < r « £ r — > = wm [ r (*■ r^) • 

where, for x,z S [0, oo), 

poo V(r) x ~ 1 zi 

T(x,z):= / u x - 1 e- u du = ^Y,— and := 0) = (x - 1)! . 

By calculating the integral, one finds that 

x~ 1 

vy ; ' j=0 J ' 

which implies that 
and that 

\ 1 / ^tan^ \ y ^ (j + y- 1)! 

P (6 < a) = 7— — jtj ( 7~~77== ) ' 



(y - 1)! \1 + y/timaj j! \ 1 + Vtan q 

Compare with the results found in [101 IT2] for x = y = 1. We also remark that, for x = L and 

y = I, U — > 1 (L — > oo), while for x = 1 and y = L, U — > — 1 (L oo). Furthermore, if x = y = L, 

then [/ — t- (L — t- oo); this follows immediately from the fact that E(T(L, p) — T(L, 1 — p)) = 0{L) 
if p 7^ 1/2, whereas the variance is also of order L. This answers a question posed in [11] (see the 
end of Section 2.2 there). 



4.2. Periodic initial configuration. We consider the following initial configuration ij for TASEP: 
choose > 1, and suppose to the right of zero we start with /c + holes, then one particle, then 

k + holes, etc. So for the configuration f) with a hole-particle pair at zero we have, for % > 1, 



r]{i) 



1 if mod(fc + + 1) = 1, 
if z'mod(& + + l) 7^1. 



To the left of zero, we have a similar, dual pattern: we start with particles, then one hole, then 
k- particles, etc. So for i < 0, 



^ ~~ \ 1 if imod(A:_ + 1)^0. 



If max{A; + , £;_} > 2, then we are in the rarefaction regime. The profile a corresponding to r] starts 
at (0,0) with a jump to the right, then k + jumps down, then a jump to the right, etc. Going to 
the left, it starts with a jump up, then k- jumps to the left, then a jump up, etc. 
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The G/M/l queue. The key to solving the distribution of 

S + :=supS^'+ and 5_ := sup S^ v _ 

n>0 ' n<-l 

is the following connection to a G/M/l queue: suppose the arrival times Ai of the queue have a 
distribution 

Ai ~ Exp{p} + . • • + Expjp} (i.i.d. k+ times) . 

The server time Bi is exponential, Bi ~ Exp{l — p}. If we define Xi = Bi — Ai, and ~ Exp{l — p} 
independent of the Xj's, then 

f 1 

S+ ~ Z + + max < , sup > JQ > . 

I J 

It is a well known fact (see for example Part II of [8]) that the second term of the right-hand side 
has the same distribution as the waiting time of the stationary G/M/l queue, which exists if 
EAt > EBi, or 

l-p 1 1 

> — p> 



p k + 1 + k + 

Note that 1/(1 + k + ) equals the fraction of particles to the right of zero. The distribution of W a 
is explicitly known in terms of the Laplace transform of Ai : define 

[s)=E(e- sAl ] 



. P + s , 

Now define A + = A + (p) as the smallest positive zero of 

z^z-<t>{{\-p){l-z)). 
Then A + is the smallest positive solution of 

(4.1) A+(l-(l-p)A + ) fc +=/+. 

Note that A+ t 1 as 1 — /fi- 1/(1 + k+). Given this constant, we know that for all s > 0, 

P(W 00 > s) = A +e - (1 - p)(1 - A+)s . 
It follows by a direct calculation of the Laplace transform that, surprisingly, 

S + ~Exp{(l-p)(l-A + )}. 
In the case of S- , we would have the arrival time 

Ai ~ Exp'{l — p} + . . . + Exp'{l — p} (i.i.d. times) , 
and the server time Bi ~ Exp{p}. This means that for 

k - 

l - p< TTkZ' 

(note that fe-/(l + is the fraction of particles to the left of zero), we define A_ = A_(p) as the 
smallest positive solution of 

(4.2) A_(l-pA_) fe -=(l-p) fe -, 
and we get that 

S- ~Exp{/>(l-A_)}. 
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This easily implies that 

(4.3) P(S_>V- 0--p)0-~ X +) 

P([7 > u) = 



/0 (l-A_) + (l-p)(l-A+) • 

Therefore 

(1 - Pu )(l - X + , u ) 
p u (l - A_ u ) + (1 - p u ){l - A +)U ) ' 
and 

p/ n < \ _ (1 - p a )(l ~ A + .g) 

1 - J p a (l-A_, a ) + (l-p«)(l -*+,«)' 
(recall (|4.1|) and (|4.2|) ). We remark that k + = k- = oo corresponds to x = y = 1 in the previous 
example. In this case, A+ = A_ = and we get (1 — p) in the right-hand side of (|4.3p . which agrees 
with the previous calculation. 



4.3. Bernoulli initial configuration. Choose p\ E [0,1) P2 E (pi,l], and assume that {f](k) : 
k > 0} is a collection of i.i.d. Bernouilli random variables of parameter pi, while {rj(k) : k < 0}, 
rj(k) is a collection of i.i.d. Bernouilli random variables of parameter pi- In this case, we see 
that the profile a performs a two-sided random walk in Z 2 : when going down-right, a starts by 
moving down, and then it has probability p\ to move down, and probability 1 — p\ to move to the 
right. When going left-up, a starts by moving up, and then it has probability p2 to move up, and 
probability \ — pi to move to the left. We can again make the connection to a queue. Using the 
notation of Section 14.21 when considering 5 + , we would get that Z + and Bi (the first step and the 
server time) are distributed as a geometric sum of i.i.d. exponentials (the number of steps between 
two right-steps, so the expected number of down-steps is 1/(1 — pi)), which is again exponential: 

Z + ,B i ~Exp{(l-pi)(l-p)}. 

Furthermore, the arrival time Ai is also a geometric sum of i.i.d. exponentials, so 

Ai ~ Exp{ptp}. 

Note that in fact, in this case we have an M/M/l queue. We need that 1 — p > p\ for S + to 
be finite. It is not hard to repeat the argument of Section 14.21 (other arguments for this classical 
problem are also possible): we find 

1 - P! 1 - p 

and from this 

S ~ Exp{l — p — p\). 

For S- we get the analogous result for 1 — p < pi by replacing p\ by 1 — p2 and p by 1 — p: 

S~ ~ Exp{p 2 - (1 - p)}- 



This means that 

P(5~ > SH 

Using that p u = (1 + u)/2, we find 

F(U>u)= . 

2(P2 - Pi) 

which is the uniform distribution on (1 — 2p2, 1 — 2pi), exactly as was found in [lOj . The reader 
can also check that the above computation leads to the same limit distribution for the angle of the 



= 1 - p - pi 
P2 ~Pi 

[l-2 Pl )-u 
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competition interface found in 



5. Computing the law of the asymptotic speed in the Hammersley process 

5.1. Periodic initial configuration. Choose A > and p > A, and assume that {v + {y) ; y > 
0) and (u~(y); y < 0) are deterministic periodic configurations of intensity A and p, i.e., u + is 
concentrated on {k/X : k > 1} while v~ is concentrated on {k/p : k < —1}. Suppose A < p < p. 
Theorem [3] then tells us that we should consider 

P ( sup{v(z) - i> p (z)} > svp{v(z) - v p {z)} ) , 

\z>0 2<0 / 

where v p is a Poisson process on R. Define 

S = sup{v(z) — v p {z)} and S~ = sup{^(z) — P p (z)}. 

2>0 z<0 

Clearly, for z > 0, we have v(z) = [A^J , where [^J denotes the integer part of x. Therefore, 

S + = sup{-v p (z) + [Xz\} 
= lsup{-v p (z) + Xz}\ 

z>0 

= -\inttv p (z) — Xz}}. 

2>0 

This infimum of a Poisson process minus a linear function is studied in |16j . Theorem 3 in [16] 
entails that for k > 

F(S + >k) = F{\mi{i?p{z) - Xz}] < -k) 

z>0 

= F(M{v p (z)-Xz}<-k) 

where p + = p + (X,p) £ (0, 1) is the positive solution of 

p+ = 1 - e- p+p/x . 

For z < we have v{z) = — \p\z\\ . So when we switch to positive z, we get 

S~ = sup{^ p (^) - [pz\} 

z>0 

= \snp{u p (z) - pz}] . 

2>0 

Now we can use Equation (7) of [16]: for k > 

P(5~ < jfe) = P(sup{P p (z) - pz} < k) 



z>0 

k 



(i-p/^)E(-^) ii V Le 



(k - i) % ^ k ^i)ip 



i 

i=0 
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In particular, P(5 _ = 0) = 1 — pj\i. Therefore, 

oo 

ns + >s-) = ^p(s-<fc) P+ (i- P+ ) fc 

k=0 

= P+ (i - p/m) E £(-P/pr ^^ fe -^(i - 

fc=0 i=0 

We were not able to simplify this formula significantly. For A = 1 and p = 2, Figure Q] gives the 
graph of P(S + > S~ ) as a function of p G (1,2). We compare it with the probability (p — p)/(p — A), 

1.0 
0.9 
0.8 
0.7 
0.0 
0.5 
0.4 
0.3 
0.2 
0.1 
0.0 

1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2.0 

Figure 1. Picture of P(S + > S~), together with uniform distribution. 

which we would get if we would take a Poisson process of intensity A for x > 0, and of intensity p 
for x < (see Section 15. 2p . 

An interesting limit is A — > oo. This means that S + = 0, and we get 

¥(S + >S-) = l-p/p. 

This is exactly the same probability as when we would have a Poisson process of intensity p left of 
0, and no particles to the right of (see Section [5T2]) ! 

5.2. Poisson initial configuration. Choose A > and p > A, and assume that {v + (y) ; y > 0) 
and {y~ (y) ; y < 0) are independent Poisson counting processes of intensity A and p, respectively. 
We will only consider the case p £ (A, p), the other cases are trivial. Define two asymmetric simple 
random walks W + and W~ , with 

F(W + (n + 1) - W + (n) = +1) = p + := 

A + p 
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and 

P 



F(W(n + 1) - W~(n) = +1) = p~ :-- 

H + p 

Since v + , Up, v~ and v~ are independent Poisson counting process, it is not hard to see that 
sup{^ + (z) — u^{z)} = sup W + (n) and that sup{z^~(z) — v~{z)} = supW^(n). 

z>0 ' n>0 z<0 ' n>0 

Furthermore, it is well known (and easy to see) that S l := sup n>0 W l (n) ~ Geo(r*), where r % 
— p 1 ), since p l < 0.5 (for i = +, — ). We find 

oo 

W(S + > S~) = J^P(S+ > fc)P(5" = k) 

k=0 

oo i _ 

fc=0 ^ 
Therefore, by Theorem [3] (and taking p„ = 1^), this proves that 



' if v > 1/A 



2 



i 



-A 



if l//i 2 < w < 1/A 2 



1 if u < l//i 2 . 
This agrees with the results found in [U [S] . 

6. Final comments 

One interesting aspect of our method is that (j2.8|) also holds in the lattice last-passage model 
with general i.i.d. weights, as soon as we can construct Busemann functions. In the exponential 
case, the limit shape is given by f(u,v) = {^fu + \/v) 2 . Also in the case of geometric weights, 
the shape function is explicitly known, and with our methods we could determine the distribution 
of the asymptotic speed of a second class particle. Busemann functions can be constructed under 
a curvature hypotheses for the respective limit shape [13] ■ However, to prove that this curvature 
hypotheses holds in a great generality is by far one the most challenging problem in lattice last- 
(or first-) passage percolation models. 

In the Hammersley last-passage model with random weights [7] (each Poissonian point x has a 
weight X x ), the limit shape (due to the invariance of the model under volume preserving maps) 
f(x, t) = j^/xi allows us to construct Busemann functions, and then prove (|2.17p for the respective 
second class particle. However, we do not know the equilibrium distribution and so we can not go 
much further than that. 

Part of the results concerning the distributional behavior of second class particles and compe- 
tition interfaces in the rarefaction regime were already known [H [91 [iOj [TTl [121 [13] . The genuine 
contributions are (|2.5p . (|2.8p and (|2.17p . and how it can be used to compute the distribution of the 
asymptotic speed, as soon as we have a good candidate for the equilibrium measure. 
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